Kryzhevich S.G. 

Saint-Petersburg State University, Russia, 
E-mail: kryzhevitz@rambler.ru, kryzhevich@hotmail.com 

Smale Horseshoe and Grazing Bifurcation in Impact Systems 



Abstract. Bifurcations of dynamical systems, described by a second order differential equa- 
tions and by an impact condition are studied. It is shown that the variation of parameters 
when the number of impacts of a periodic solution increases, leads to the occurrence of a 
hyperbolic chaotic invariant set. 

Introduction. The vibro-impact systems appear in different mechanical problems (mod- 
eling of impact dampers, clock mechanisms, immersion of constructions, etc.). All the impact 
systems are strongly nonlinear. Their properties resemble a lot ones of classical nonlinear 
systems. Particularly, the chaotic dynamics is possible [1] - [5], [7] - [9], [12] - [15], [19], [20]. 

There is a big number of publications, devoted to bifurcations, proper to vibro-impact 
systems. One of them, the so-called grazing bifurcation, first described by A. Nordmark 
[15], corresponds to the case, when there is a family of periodic solution, which has a finite 
number of impacts over the period and this number increases or decreases as the parameter 
changes. For the bifurcation value of the parameter, the periodic solution has an impact with 
a zero normal velocity. It was shown that this bifurcation implies the nonsmooth behavior 
of solutions, instability of the periodic solution in the parametric neighborhood of grazing 
and, in additional assumptions, the chaotic dynamics (see [4], [5], [7]- [10], [15] and [20] 
and the references therein). However, there was no analytical conditions, sufficient for the 
existence of chaotic invariant sets. For the s.d.f. case these conditions have been obtained in 
the author's work [13]. In this paper generalize this result for systems with several degrees 
of freedom. The main result of this paper is the method, which allows to find homoclinic 
points, corresponding to grazing. The main idea of the proof is the nonsmoothness of Perron 
surfaces in the neighborhood of periodic solution. If these manifolds bend in a "good" way 
(the corresponding sufficient conditions can be written down explicitly) they can intersect. 
This implies chaos. We study a motion of a point mass, described by system of second order 
differential equations of the general form and impact conditions of Newtonian type. In order 
to avoid technical troubles we assume that the delimiter is plain, immobile and slippery. 
However, there is no obstacles to apply the offered method to the systems with a mobile 
delimiter (see, for example [9]), ones with non- Newtonian model of impacts [2], [7], [10], [11] 
and even to some special cases of strongly nonlinear dynamical systems without any impact 
conditions. 

1. Mathematical model. Consider a segment J = [0,/i*] and a C 2 smooth function 
f(t,z,/i) : M 2n+1 x J 4 1" Suppose that f(t,z,n) = f(t + T(/i), z, //). Here the period 
T(fi) is a C 2 smooth function of the parameter /i, and T(0) > 0. We may suppose without 
loss of generality that T(/z) does not depend on /z, making, if necessary, the transformation 
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f = tT(0)/T(n). Denote 
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Consider the system of In first order ordinary differential equations 

x k = Vk] Vk = fk(t,z,fj,), fc = l,...,n. (i) 

Let r = r(yi, [/,) G (0, 1] be a C 2 - smooth function. Let us note by col(ai, . . . , a m ) the 
column vector, consisting of elements Suppose that the system (1) is defined for 

z G A = [0, +oo) x M 2n_1 and the following Newtonian impact conditions take place as soon 
as the first component x 1 of the solution vanishes. 

Condition 1. 

1. If Xl (t ) = then x(t + 0) = x(t - 0), yi (t + 0) = -r(- yi (t - 0),fi)yi(t - 0), 
y(t + 0) = y(t -0). 

2. Let the solution z(t) = col(zi(t), . . . ,z n (t)) of the system (1) be such that zi(to) = 
for a certain instant to and /i(io, 0, z(to),fj) ^ 0. Let ((t) be the solution of the system 

x k = yk, y k = f k (t,0,0,z,fj,), k = 2,...,n 

with the initial data ((t ) = z(t ) and the segment I = [t ,ti] be such that 

/i(*,0,0,C(*),A*)<0 

for all t G /. Then Zi(t) = and z(t) = ((t) for any t G i\ 

Remark 1. Similarly one may consider a domain A = x IR n , where the domain A x is 
bounded by a C 2 smooth manifold M of the dimension n — 1. 

Let us denote the dynamical system, defined by Eq. (1) and mentioned impact conditions 
by (A). 

2. The family of periodic solutions. Since the solutions of the vibro-impact systems 
are discontinuous at the impact instants, the classical results on integral continuity are not 
applicable. Nevertheless, the following two statements hold true. 

Lemma 1. Let z(t) = col(xi(t), yi(t), . . . ,x n (t),y n (t)) be the solution of the system (A) for 
fj, — (j, and the initial data z(t ) = z° = col(x1, . . . , y°) and x\ ^ 0. Suppose that 
this solution is defined on the segment 3 t . Assume that there are exactly N zeros 

t- < rf < . . . < r^r < t+ of the function x\(t) over the segment [t-,t+] and yi(Tj — 0) ^ 0, 
(j — 1, . . . , N). Then for any e > there exists a neighborhood U of the point col (to, z°, Ho) 

N 

such that for any fixed t G [£-,£+] \ U (t° — e, rj} + e) the mapping z(t, t±, z 1 , fj,) is C 2 - 

k=i 

smooth with respect to the variables (t^z 1 ,//) G U. Moreover, this solutions have exactly N 
impact instants Tj(t\, z 1 , fi±) (j = 1,...,N) over the segment [t_,t + ]. These instants and 
corresponding velocities Yj = — yi(Tj(ti, z 1 , /ii) — 0, t±, z 1 , fj,i) C 2 smoothly depend on t x ,z x 
and Hi. 



Proof. Let the number k be such that t e [rk,r k+1 ] (assume, if necessary r = i_, 
Tjv+i = The solution of the vibroimpact system (A) is also one of (1) over the 
segment [r fc ,r fc+1 ]. The impact instants r k and r fe+1 as well as the impact points 

z(r k - 0,£i,z 1 ,//i) and z(T fe+ i - 0, ii, z 1 , /ii) 

smoothly depend on their parameters. Similarly, the values r k -i and z(rfc_i — 0, ti, z 1 , Lii) 
smoothly depend on r k and z(r k — 0, ti, z 1 , Lii), as well as Tfc +2 and z(r k +2 — 0, £i, z 1 , /ii) are 
C 1 - smooth functions of r k+1 and z(r k+ i — 0, t 1 ,z 1 ,/ii) and so on. This proves the Lemma 
1. 

Condition 2. (jf£</. ij. T/iere exists a continuous family ofT - periodic solutions 
(p(t,fi) = col{ V x {t,Li),p y {t,Li),...,^ x n {t,Li),^ y n {t,Li)) tet, i*e J 
of the system (A), satisfying the following properties. 

1. For all ji ^ the component <p x (t, //) /ias exactly N + 1 zeros r (/i), . . . , Tjv(/i) over £/ie 
period [0, T). 

T/ie velocities Y k ([i) = —(p\(r k (fi) — 0,/i)) are stic/i t/iat 

Y^aO > /or a// /i > 0, F o (0) = 0, 

A(r o (O),O,O,v9(r o (O),O),O) = 0o>O, (2) 
Y k (fi)>0, VlieJ, k = l,...,N. 

Here <p(t,n) = coh>| (t, Li), ip y 2 (t, //),..., ip x (t, Li), ip y n (t, li)). 




3. The instants r fc (/i) and the velocities Y k (fj,) continuously depend on fi e J. 
Denote 

ifi X (t, Li) = C0l(^f(t, Li), . . . , <p$(t, Li)), if y (t, fi) = COl(^?(f , Li),..., tp V n(t, Li)), 

<f x (t, li) = col((p x (t, /i), . . . , (f%(t, /i)), <p»(t, li) = col(^(t, Li), . . . , ip y n (t, fx)). 
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We may suppose without loss of generality that r (/i) = 0, </?(0, 0) =0. This may be obtained 
by the transformation t' — t — r (/x). Define 

9 = — min (min Ti(/i),T — maxT^/i) 

Fix a small positive 9 < 9 and consider the shift mapping for the system (A), given by the 
formula S^g(z°) = z(T — 9 + 0, —9,z°,ii). For small positive \i and 6 the mapping S^g is 
C 1 - smooth in a neighborhood of the point z^g = ip(—9,ii). Without loss of generality, we 
may assume that 

lim Za,e = 0. 

3. The separatrix. Denote T^g = {z° G A : 3ti G [-T,T] : z^h, -6,z°,n) = 0}. 

Lemma 2. There exists a neighborhood Uq of zero such that if the parameters fi and 9 are 
small enough, the set r Mi # Q Uq is a surface of the dimension 2n — 1, which is the graph of 
the C 2 smooth function x\ = r y fJ , t g(x,y) (fig. 2). Moreover, 



7^M = yn /lM ; w) +7,,«M), (3) 



where 7^0 is a C 2 smooth function such that 7 M ,e(0) = 0. 




Proof. Take a point ( G T^g. Let the instant to be such that Zi(t , —9, () = 0, s = t — t , 
Zi(t + 0, — 9,() = col(xi(t),yi(t)). Let us show that if to is close enough to —9, we may 
take s ^ l^o + 9\ so that the function Xi(t + s) does not have zeros on [— s ,s ], except 
s = 0. Otherwise, there exists a sequence t k — > —9 (suppose without loss of generality, 
that tg > —9 and the sequence decreases), a sequence t\ G [—0,$) and one, consisting of 
solutions, uniformly bounded on the segment [— 0, maxi^]: 

z\t) = col(^(t), . . . , 4(t)) = co\(x k (t), y*(f), x k n (t),y k n (t)) 
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of the system (A), such that = 0, x\{t\) = 0. Also, there exist time instants 



t\ G (ti,to), that xXit*) = and instants £3 G (t^o) sucn that ^1(^3) = 0- Moreover, 
t\ — > —6, Xi(tg) — > 0, — >■ 0. Without loss of generality, we assume that z(t§) — > z°. 

Then £i(tf) ->■ ^(-6,2°,^) = 0. This contradicts to (2). 

Then for all s G [— s , s ] the function Xi(t + s) can be presented as series 

si (to + s) = X 2 s 2 + X 3 s 3 + . . . (4) 

Differentiating (4), we obtain that xi(to + s) — 2X 2 s + 3X^s 2 + . . .. On the other hand, 
X 2 = %(t + 0)/2 h{-Q,z^^)l1 as t -0- Then 

= /i(-^^,/i)(to + ^) 2 (l + o(l))/2; y 1 (-0) = / 1 (-0,z Mi< ,,/i)(* o + 0)(l + o(l)). 

Since yi = £1, the formula (3) is true. The lemma is proved. 

Take a small parameter q > such that the sets V^g = {z G A : ||z — z^ e \\ ^ ^} C £/ , 

^ M 7e = {(&> 2/) e v nfi '■ x i < Infill y))i V £e = {(x,y) e V : X! > 7^0(2, y)} 

are correctly defined and nonempty. 

4. The main result. Consider the matrix 

oz 

A= lim — (T- 6 + 0,6, z°,Q) 
9^0+ dz° 

Let Ao = det A. Denote the elements of the matrix A by ay and ones of the matrix A~ x by 
ay. Denote the columns of matrices A and A 2 by Aj and A 2 respectively, the strings of the 
matrix A^ 1 by Aj. If n > 1 and ai 2 7^ consider the (2n — 2) x (2n — 2) matrix A = (ay), 
defined by formulae 

dy = h a i+2 j+2, = 1, . . . ,2n - 2). 

01a 

Similarly, if n = 1 and «i2 7^ 0, we define the matrix A = (ay) 

ay = h a i+ 2j +2 , («, J = 1, • • • , 2n - 2). 

ai2 

Assume that the at least one of the following statements is true. 
Condition 3. 



1. Either n = 1 or the matrix A does not have eigenvalues on the unit circle in C. 
2. 

2n 

ai2 > 0, ^2a lk a k2 < (5) 
k=i 

Condition 4. 

1. Either n = 1 or the matrix A does not have eigenvalues on the unit circle in C. 
2. 

2n 

«i2 > 0, ^2a lk a k2 < 0. (6) 

k=l 

From the geometrical point of view, the first items of Conditions 3 and 4 mean that 
the fixed points z^g are saddle hyperbolic and the inequalities (5) and (6) provide that the 
corresponding stable and unstable manifolds intersect. This will be shown below. 

Later on we shall suppose that Condition 3 is satisfied. Otherwise, we consider the map- 
ping instead of S^ t g. Then the matrix A is replaced with A -1 , and the condition (6) 
with (5). The similar reasonings shall prove the statement of the lemma in the considered 
case (see the right part of the figure 4). All the proofs given below, may be repeated for this 
case. 




Theorem 1. // Condition 2 and one of Conditions 3 or 4 are satisfied, there exist values 
fiQ > and 9 > such that for all fi G (0,/io) the mapping S^o is chaotic in the sense of [6J. 
More precisely, there exists an integer m and a compact set K = K^ e invariant with respect 
to S™ 9 such that the following conditions are satisfied. 

1) There exists a neighborhood U^ t g of the set K^g such that the mapping S™ d \u^ e is a 
diffeomorphism. The invariant set K„ 9 is hyperbolic. 
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2) The mapping S™ d \K^ e has infinitely many periodic points. 

3) The periodic points of S™ 6 are dense in K^g. 

4) The set K^g is transitive i.e. there exists a point p^g e K^g, such that the orbit 

{sfyipM^-oo ls dense m K »fi- 

Remark 2. The similar results may be obtained for the systems with any finite number of 
grazings over the period. 

5. Grazing. Now we start to prove the theorem. Note that all the mappings S^g, corre- 
sponding to the same value of //, are conjugated. Fix a small value \i > and a solution 

z°(t) = C o\(x l (t)y 1 (t),...,x° n (t)Mt)) 

of the corresponding system, having an impact at the instant to- Suppose that the corre- 
sponding normal velocity Y 0i = —y®(t — 0) is nonzero. We consider Yn as a small parameter. 
Denote Z = z°(to — 0). Fix a positive value so and consider the mapping 

G(C) = z(t + s , *o - s , C, a0, 

defined in a neighborhood of the point (° = z (t — so). Here we assume that the point (° and 
the parameter s are chosen so that there exists a neighborhood Q 3 (° such that any solution 
z(t) = z(t, t — s , Z-,fj) (z- G Q) has exactly one impact over the segment [t — s , t + s ]. 
Denote the corresponding instant by t± = ti(z-) and the normal velocity, defined similarly 
to Yoi, by Y\ = Yi(z-). Let Z = z{t\) be the tangent component of the solution z(t) at the 
impact instant. Take the values s± = s±(z_) so that to±s = Ai(z-) ±s±(z_) for all z_ e Q. 
The mapping G is smooth in the neighborhood of the point Co, let us estimate the Jacobi 
matrix DG(( ). Denote 

z + = z(to + s ) = z(t + s , h + 0, 0, rY 1 , Z, p), 
x + = x(t + s ) = x(t + s , h + 0, 0, rY ± , Z, fj), 
V+ = z(t + s ) = y(to + s , t x + 0, 0, rY 1} Z,ji), 
x- = x(t - s ) = x(t - s , h - 0, 0, -Y u Z, p), 
y- =y(to-so) =j/(i o -so,*i-0,0,-Yi,Z,^). 

Similarly, we define the values x lt ±, yi t ±, x±, y±. Consider the Taylor formula for values 
z± as functions of s±: 

Xl - = Fis_ + fi(t u 0, -Y U Z, ii)s 2 _/2 + p xl 4s-,h, Fi, Z, p)sl; 
Vi- = ~Y\ - /i(*i,0, -Yi, Z, fi)s- + p yl -(s_,t 1 ,Y 1 , Z , n)s 2 _] 

x. = x(h) - y(h - 0)s_ + f(h, 0, -Y u Z, /i)s 2 „/2 + p s __(s_, h, Y u Z, /i)s 3 _; 
V- = V(ti - 0) - f(t 1 ,0,-Y 1 ,_Z,fi)s- + py-(s-,t 1 ,Y 1 ,Z,ji)sl.; 

xi, + = rYiS + + /i(ti, 0, rYi, Z, p)s\/2 + p x i+(s+, t u Y u Z, p)s 3 + ; 



(7) 



Vi,+ = rYi + fi(t 1 ,0,rY 1 ,Z,n)s+ + p yl+ (s + ,ti,Yi,Z,fi)s 



2 

+ ' 



x + = x(h) + y{h - 0)s + + /(_*!, 0, rF 1? Z, ^)4/2 + p s +(s+, *i, Fi, p,)s 3 + ; 
V+ = V(ti - 0) + /(tx, 0, rF 1? Z, /i)s + + p s+ (s + , t 1? F 1? Z, /i)s^. 

Here all functions, denoted by the letter p with different indices, are C 2 smooth with respect 
to all arguments except s±. Denote 

- - „ <9r 

fok+ = fk(to, 0, rloi, Z , fj), f ok _ = f k (t o ,0, —Y 01 , Z , [j,), r = r + — — F i- 

°yi yi =Y i 
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It follows from (7) that 



dz. 



+ 



d(s + ,Y l7 Z) 



d(s-,Y 1: Z) 



s+=0,Y 1 =Y i,Z=Zo 



s-=O,Y 1 =Y 01 ,Z=Z Q 
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E 2n -2 



Here E 2n -2 is the unit matrix of the corresponding size, 

Q+ = co\(y 2 (t - 0), / 02 +, • • • , y n (t - 0), f 0n+ ), 
Q- = co\(-y 2 (t - 0), -/02+, • • • , -y n (to ~ 0), ~fon+)- 



Denote 



f'k — 



df k (t ,0,y 1 ,Z ,fx) 



dyi 



2/1=0 



— r 









&21 


— r 





(8) 


#1 





E2n-2 ) 





Clearly, ds + /ds_ = — 1. Then, similarly to the results of the paper [10], we obtain 

B = Km ^± = 

(fig. 2). Here £?i = col(6 3 i, . . . b 2n i), 

hi = -(/oi+ + r/oi-)/loi = -(r + 1)0 O (1 + O(Y 01 ))/Y 01 , b 2j -n = 0, 
b 2j i = (fo j+ ~ fa-)/Y 0l = (r + l)f k + O(Foi). 

Note that, det B = rr — >■ r 2 as F i — >■ 0. 

It is shown on the figure 2 how the small neighborhood R of the form 

R = {(6, Vi, • • • , £n, VV.) : & " C?l ^ 5, 1^ - < 5} 
is stretching under the action of the mapping Si^g, defined by the formula 

s 1>lifi (0 = z(e + o,-e,{,ri. 

Here R 1 = S 1>fJ , >0 (R). 

6. Lyapunov exponents. We check that the fixed point z^^ of the mapping S^q is hy- 
perbolic and estimate the bigger and the smaller absolute value of eigenvalues of the matrix 
D = DSp^Zpfi) and ones of small perturbations of this matrix. For a fixed value \i > denote 
D = DS^(z^). The mapping S^q can be presented as the composition S^g — £2,^,0 £1,^,0, 
where S 2 ^e(C) = Z (T — 9 + 0, 9, (, [/). The Jacobi matrix A^ e = DS 2 ^fi(Si^fi(z^ e )) tends to 
A as fi, 9 — > 0. The matrix = DSi^q(z^q) is of the form (8), where Y 01 — (/■*)■ Then 

L> = A^B^ = (-(r+l)A 2 <f> (l+O(Y ))/Y , -rA 2 (l+O(Y )), A 3 +O(Y ), A 2n +O(Y )). 



Since det D = det A^ g det f?^ 



0(1q))Ao, if the inequalities (5) take place and if fi 



is small enough, one of the eigenvalues of the matrix D is 

A + = -(r + l)a 12 0o(l + O(^o))Ao- 



The corresponding eigenvector u + equals to A 2 +O(Y ). The eigenvalue A + is of the multiplic- 
ity 1, the linear space, corresponding to other eigenvalues, tends as jj, — > to the hyperplane 
7Ti, given by the condition X\ = 0. Since a\ 2 ^ 0, the vector A 2 is out of 7Ti. The matrix D^ 1 
satisfies the following asymptotic estimate 

/ -rA 1 + O(Y ) \ 



(r + l)Ai(l + O(Y Q ))/Y 
A 3 + O(Y ) 



\ 



J 



A 2n + O(Y ) 

It follows from the form of this matrix, that the matrix D~ x has the eigenvalue 

A: 1 = (r + l)0 o ai 2 (l + O(Y ))/(r 2 Y ). 

The corresponding eigenvector satisfies the asymptotical estimate w_ = e 2 + 0(Yq). 

If n — 1, it is clear that the matrix AB as well as matrices DS^gi^z) corresponding to 
points z of a small neighborhood of z^g are hyperbolic. Otherwise we need the following 
lemma. 

Lemma 3. If Condition 3 is satisfied there exist positive constants /i and eo such that for 
any matrix A' such that \\A' — A\\ < e and any /i G (0,// ) the matrix A'B^g does not have 
eigenvalues on the unit circle in C. 

Proof. Suppose the statement of the lemma is not true. Then there exist a sequence of 
matrices Ak —> A and sequences /ik, 8k —> such that all matrices A k B fJikt g k have eigenvalues 
of the form \ k = exp(iip k ). Denote the corresponding eigenvectors by u k = co\(u k , . . . , u 2n ). 
Without loss of generality we may assume that \\u k \\ = 1 for all k and that 



Afc ->■ A = exp(i^o) 



u k ^u° 



col(w°, 



> U 2n) 



as k — > oo. Denote the elements of the matrices Ak and B fJikt g k by a k j and respectively 
(A; e N, i,j = 1, . . . ,2n). Let r fc = r(Y (/j, k ), fj, k ), h = r(Y (fi k ), 

2n 



Note that since a^ 2 > 012/2 > 0, b\ 2 — > 00 and other values b\ 2 are uniformly bounded, the 
sequence d\ tends to +00 as k — > +00 and 

lim d k jd\ = a i2 /a 12 , (i = 1, . . . , 2n). 

k— >+oo 

From the definition of eigenvalues and eigenvectors we obtain 



2n 



d k u k - r k a k 2 u k + a js u s = A*uJ j = l,...,2n 



(9) 



s=3 
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It follows from the first equation (9) that 

2n 

-r k a 12 u\ + Yl a is u s 

-d\ + X k ^° 

as k — > +00. Consequently, u\ = 0. Substituting this expression to Eq. (9), corresponding to 
j > 1, we have 

2n 

d){-r k a\ 2 u\ + £ a^u*) 2 « 

s=3 r fc aj 2 u£ + ^ aj s uj = A fc ^ fc j = 2, . . . , In. 

til + A fe s=3 

Proceeding to the limit as A; — )■ 00 we obtain 

2n 



^2 + a,. U; = A„»», j = 2,...,2„. 

&12 



Clearly, at least one of values u® (j — 3, . . . , 2n) is nonzero. Then A is an eigenvalue of the 
matrix A. This contradicts to our assumptions. 

7. Homoclinic point. The mapping S^g is differentiable at the points of the set V~ . Due 
to Lemma 3 the eigenvalues of Jacobi matrices D are out of the unit circle, provided /i is 
small. Then, due to the Perron theorem, in a small neighborhood of the point z^g there exist 
the local stable manifold W s and the unstable one W u of the mapping S^g. Both of them 
are smooth surfaces. Let M s ' u = T z ^ e W s ' u , n s = dimW^, n u = dimW u = 2n — n s . Select 
the orthonormal basises e\ , . . . , e s n and e", . . . , in the spaces M u and M s respectively 
so that e\ = u-/\\u-\\ and e" = tx + /||w + ||. Extend the stable and unstable manifolds up 
to the invariant sets W s and W u of the mapping S^g. The obtained sets consist, generally 
speaking, of a countable number of the connected components. Every of these components 
is a partially smooth manifold. 

Lemma 4. The manifolds W s and W u intersect transversally at a point p 7^ z^g (fig. 4)- 

Proof. If conditions (5) be satisfied then, for small values of the parameter ji the man- 
ifold W u intersects transversally the surface T^ e f] V^g. Denote the manifold, obtained in 
the intersection, by q u . The dimension of q u equals to n u — 1. Denote = S^g(q u ). The 
neighborhood V^g of the point z^g may be chosen so that both the manifolds q u and 
intersect with V^g and the intersection W u f] V^g \ q u consists of two connected components. 
Denote one, which does not contain the point z^g, by L u . The neighborhood V^g may be 
chosen so that diamV^ — > as /i, 9 — > 0+. Let L\ = S^g(L u ). For any z G L u the tangent 
space M u (z) = T Z W U is the linear hull of unit orthogonal vectors E%(z), E%(z), . . . , E% u (z), 
which can be chosen so that for any k — 1, . . . , n u 

lim max \\Et{z) - eV\ = 0. 

/i,0->O+ zeL u 

The surface W u is not smooth in the neighborhood of the manifold q™. For the points z G L\ 
the tangent space M u (z) = T Z W U is the linear hull of unit vectors E™(z), E% (z), . . . , E% u (z), 
where 

lim max \\E%(z) - et\\ =0 
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for all k > 1 and 

lim max\\E?(z) - AZ/\\AZ\\\\ = 0. 

ti,e^o+ z^li 1 2/11 211 

It follows from the conditions (5), that for any z G L\ the set E\(z), E 2 {z), . . . , E™ (z) 
is linearly independent. The vectors A 2 and A\ lie at different half-spaces, separated by 
the hyperplane ~K\. Consequently, for all z G L u , Z\ G L\ and z 2 G W u f] V^g the vec- 
tors Ef(z ) and E™(zi) lie in different half-spaces, separated by the linear hull of vectors 
E% (z 2 ), ■ ■ ■ , E™ (z 2 ), E((z 2 ), . . . , (z 2 ). This means that for the small values of the param- 
eter fj, and 9 the neighborhood V^q may be chosen so that the surfaces L\ and W s intersect 
transversally. This proves the lemma for the considered case (see the left part of the figure 
4). 

The Smale-Birkhof theorem [18] on the existence of a chaotic invariant set in a neigh- 
borhood of a homoclinic point is not applicable in the considered case since the mapping 
is discontinuous. However, the similar techniques will help us to find a chaotic set of the 
mapping S^e- 
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ABSTRACT 
Kryzhevich S. G. 

Smale horseshoes and grazing bifurcation in vibroimpact systems 

Bifurcations of dynamical systems, described by a second order differential equations and 
by an impact condition are studied. It is shown that the variation of parameters when the 
number of impacts of a periodic solution increases, leads to the occurrence of a hyperbolic 
chaotic invariant set. 



